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Abstract 

A finite difference computer model is developed to simulate the exposure statistics of a radio 
freciuency buoyant antenna as it is towed in a random seaway. The model allows th(^ user 
to iH'escribe antenna properties (length, diameter, density, etc.), sea conditions (significant 
wave height, develoi)inent of sea), and tow speed. The model then simulates the antenna-sea 
interaction for the desired duration to collect statistics relating to antenna performance. The 
model provides design engineers with a tool to predict antenna performance trends, and conduct 
design tradeoff studies. The antenna envisioned is a submarine floating antenna which would 
(uiable coinmunications at speed and depth, greatly enhancing the stealth and survivability of 
the US Navy’s submarine force. 
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Chapter 1 



Introduction 



Driven largely by deiriaiid in the corninercial sector, the capability for high speed data transfer 
has grown at an amazing rate over the last few years. Recently deployed military communi- 
cations systems have been able to leverage off their commercial counterparts to some extent, 
and now operate at higher frequencies and data rates than the systems they have replaced. 
However, the United States submarine force faces unique challenges in fully utilizing systems 
which operate at higher frequencies. Submarine operations are typically driven by the need 
for stealth, and stealth is best preserved by remaining fully submerged. Currently submarine 
communications can be typified by a pattern of coming to periscope depth periodically to raise 
a radio mast, send or receive required message traffic, then submerge. While submerged, radio 
comniunications systems are limited to low data rate receive-oiily modes transmitted at very 
low frequencies (ELF or VLF bands). Clearly, the communications goal of the submarine force 
is to be able to utilize the high frequency (UHF, SHF, or EHF bands) and high data rate 
conimunication circuits while retaining the stealth of operation at depth. 

This goal of high data rate communication at depth was the motivation of a U. S. Navy Ad- 
vanced Technology Demonstration (ATD). Initially funded jointly by DARPA and ONR in FY 
96-97 as a feasibility study, the Buoyant Cable Array Antenna research became a DARPA/ONR 
Demonstration in FY 98. The program was promoted to the level of Navy ATD in FY 00. This 
ATD. one of only 6 ATD’s funded by the Navy for fiscal year 2000, reflects the high priority 
the Navy places on submarine connectivity. The goals of the program include [1] 
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• verify the capability to conduct beamforining signal processing 

• verify adequate antenna perforiiiaiice and gain in a seaway 

• develoj) traiisiiiit capability 

• conduct design studies on key antenna characteristics (diameter, length, density, etc.) 

Ill addition to the engineering challenges inherent in any new design work, the program 
managcn- must ensure the designers keep system affordability and flexibility in mind. 

1.1 Background 

Lincoln Laboratory is a federally funded laboratory administered by Massachusetts Institute 
of Technology. Lincoln Lab, specializing in radar and electrical engineering, is well qualified to 
tackle the challenge of designing a floating linear array RF communication antenna. However, to 
better characterize the behavior of the antenna in a random seaway, Lincoln Lab added Professor 
Jerome Milgram of the MIT Ocean Engineering department to the antenna development effort. 
Professor Milgram was made the head of the Exposure Statistics Working Group, whose charter 
was to develop a model to provide statistics on 

• the average fraction of antenna elements exposed over an extended period of time 

• the temporal variation of the exposure statistics, including duration of reception dropout 
times 

• the spatial distribution of those exp)osed elements 

This thesis work was undertaken to support the charter of the Exposure Statistics Working 
Groiij). By constructing a finite difference comi)uter model of the antenna, a design tool was 
envisioned that would allow engineers to evaluate performance trends and pose “what if'* design 
(piestions. And most importantly, the computer model would permit performance trade-off 
studies to be conducted at the system engineering level. 
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Chapter 2 



Model Development 



2.1 Characterization of Sea 

2.1.1 Linear Theory Ocean Waves 

For the j)urposes of this work, the sea was assumed to be described by linear theory, random, 
deep-water waves. “Deep water” is defined to exist when X < 2d where d is the water depth. 
The random aspect of the waves implies that 

• the distribution of wave elevations from the calm water level follow a gaiissian form, and 
hav(^ zero mean 

• the sea statistics obtained from a single location are the same as the statistics obtained 
from an ensemble of locations (ergodic) 

• the sea statistics remain constant with respect to time (steady state process). 

Because the last condition is so restrictive, we relax it to allow the statistics to vary slowly 
with time (weakly steady state process). A summary of the deep water linear theory wave 
relations is provided in Table 2.1. 

Assuming a lower sea spectral energy bound of 0.7 rad/sec, which is reasonable for the 
range of seas analyzed in this report, this translates into a deep water depth of 252 meters. 
For depths shallower than about 250 meters, the deep water theory begins to break down, and 
a de[)th-dependent treatment of the sea waves must be considered. Referring to the formulas 
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Quantity 



Formula 
k — 



Wave Number 
Circular Frequency 
Dispersion Relation 
Velocity Potential 
Surface Displacement 
Radian Frequency 
Phase Velocity 
Group Velocity 
Horizontal Fluid Velocity 
Vertical Fluid Velocity 
Wave Energy Density 



k=^ 

0 = siii{kx' — Lot) 
h = acos(fcx' — Lji) 

u = y/^- 

Cg = 5 y/^' 

V = au>e^' cos{kx' — u>t) 
tv = awe^' sin{kx' — uit) 

Z7* 1 2 

E = ^pga- 



Note: .x' refers to a stationary x coordinate system. 



Table 2.1; Linear Theory Deep Water Wave Relations 



ill Table 2.1. the wave forms in linear theory can be seen to be sinusoidal. However, simple 
observation of the sea surface shows that the wave forms display a complicated structure which 
has random features. This observed behavior can be approximated by adding together a series of 
sinusoids, each frequency component having the appropriate amplitude and a random phase. As 
the number of frequency components increases, the approximation approaches the continuous 
distribution of the actual sea. This discrete method of describing a one dimensional random 
seaway is presented in equation 2.1 



N 

h{x\t) = ^a;(c(;)cos(A’ 2 (a;)x' (2.1) 

t=i 

where 

h{x\t) is the surface elevation of the sea, referenced to the calm water level 
a^(u;) is the wave amplitude of each discrete frequency component 
kt{uj) is the wave number of each frequency component 
is the set of discrete circular frequencies chosen 
is the random phase assigned to each frequency component 
and the term has been dropped, since evaluation occurs at the sea surface {z = 0). 

It is fairly straightforward to determine most of the terms on the right hand side of equation 
2.1: k'i is deteriiiined from the dispersion relation, (f>^ is a random number with a uiiiform 
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Figure 2-1: Sample Ocean Wave Energy Density Spectrum: Hs = 2.0 m 

distribution [0, 27 t]. The evaluation of the amplitudes, however, requires some discussion. 

We begin by defining the spectral energy density function, 5(a;), which represents the energy 
of the ocean waves having frequency u. (See Figure 2-1 for a typical energy density spectrum.) 
The units of S{u)) are defined to be in vi^ / sec, so that when S{u) is integrated with resj)ect 
to a;, the result is energy in units of m^. (All that is lacking to give this quantity the normal 
units of energy is the ^pg term seen in Table 2.1 for the wave energy density.) The spectrum of 
Figure 2-1 is said to be single-sided, which refers to the fact that the frequencies are constrained 
to l)e positive. For a double-sided spectrum (positive and negative frequencies allowed, and the 
spectrum defined such that S{u)) = S{—u))) of equal sea severity, one must divide the spectrum 
by 2, since for constant sea severity the area under the spectral energy density curve must be 
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Figure 2-2: Determination of Component Amplitudes cij{u;) 
constant. The amplitudes a^(cj) are obtained from Eq. 2.2 

= 5(cj)Aw (2.2) 



where 

(u{io) is the amplitude of the band of sinusoid waves Auj with center frequency uj. 

Eq. 2.2 is really nothing more than a statement that the square of the wave amplitude is a 
measure of the average energy of the waves, because both sides of the equation represent energy 
in units of nr. A method of obtaining the amplitudes ai{u) is shown in Figure 2-2. 

Although S{lj) has been presented in terms of wave energy content at a given frequency, it 
can be shown from the Weiner-Kintchine theorem that S{u) is related to the wave amplitude 
record at a fixed location in the ocean. Following Ochi, [2] we define a wave record at a given 
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location as Then the autocorrelation of x'{t) is given by 

R{r) = Vim x'{t)x'{t + r)dt (2.3) 

7 — oo J-T 

and by the Weiiier-Kint chine theorem, the spectral energy density function, 5(cj), and the 
autocorrelation function. i?(r). form a Fourier transform pair. 



5(a;) = 


1 


roo 

/ R{t) 


(2.4) 


7T „ 


I — OO 






roo 




R{t) = 


1 5(a») 

— OO 


(2.5) 



This method’s value lies in its ability to convert observed wave records into spectral energy 
density functions. 

2.1.2 Spectral Energy Density Functional Forms Assumed for the Model 

Numerous forms have been proposed by researchers attempting to characterize the sea wave 
environment [3] [4] [5] [6] [7]. The model incorporates the spectral formulation as a subroutine, 
hence it is completely arbitrary from the model viewpoint how the spectrum is defined. In order 
to limit the scope of the simulation parameters, yet retain the ability to vary the sea energy 
form, two wave spectral energy subroutines were developed for the model. Though they may 
appear very different at first glance, most spectral formulations share the basic framework of 
Eq. 2.G 

S{u>) = (2.G) 

wliere 

A, D are constants to be determined 

Hs is the significant wave height, which is defined to be the average wave height of the one 
third highest waves. Sea severity is most commonly given in terms of significant wave height 
among researchers. However operators (such as U.S. Navy personnel) are frequently more 
familiar with sea state. (Appendix E provides a conversion between the two systems.) The two 
spectral formulations chosen for the model are the Bretschneider two parameter spectrum, and 
the Ochi six parameter spectrum, given in Eqs. 2.7 and 2.8 respectively [2]. 
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(2.7) 




Bretschneider Two Parameter Spectrum 

The two parameters of the Bretschneider spectrum (Eq. 2.7) are significant wave height. 7/^. and 
modal frequency, Um. Modal frequency is that frequency at which S{u) reaches its maximum 
valiK'. and significant wave height is as previously defined. For the special case of a fully 
developed sea (also referred to as a fully arisen sea), i.e. one in which the wind imparts energy 
to th(‘ waves at a rate equal to that at which viscous damping dissipates their energ}' and waw- 
wave interactions transmit energy from one frequency to another; a relationship exists between 
Urr^ and Hg. Designating this modal frequency as have 

= 0.4yj (2.9) 

which establishes the energy equilibrium characteristic of a fully arisen sea. The two require- 
ments for a fully arisen sea are 

• sufficient fetch (i.e. the distance over which the wind is blowing over the water.) 

• sufficient duration (i.e. the length of time the wind has been blowing.) 

If the fetch or duration are insufficient, the modal frequency will be greater than that given 
by Ecj. 2.9; if the winds are subsiding, the modal frequency will be less than that given by 
E(j. 2.9. This is best understood by considering the mechanism by which waves are generated. 
The wind causes ripples (capillary waves) on the sea surface, which is a short wavelength, high 
frequency phenomenon. Thus as a storm is building, wave energy input occurs at high freciuency, 
and the energy gradually is transferred through wave-wave interactions to lower frequencies. 
Likewise, as a storm is subsiding, the high frequency waves attenuate most rapidly due to 
viscous effects, and the result is a skewing of the spectral energy toward lower frequencies. By 
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Parameter 


Value 


Hsi 


OMHs 


Hs2 


O.MH, 


^’ml 




^7712 


1 15p— 0.039Hs 


Ai 


3.00 


A2 


1.54e-o.062W, 



Table 2.2; Ochi Six Parameter Spectrum: Most Probable Values in Terms of Hs 

introducing a term we will call the development Dev. we can quantify the degree of sea arousal 
as given in Eq. 2.10. 

~ Dev • ( 2 . 10 ) 

Using this convention, values of Dev <1.0 correspond to a decaying sea. and values of Dev > 
1.0 correspond to a building sea. Figures 2-3 and 2-4 show representative plots of S{lv) as 
developinent and significant wave height, respectively, are varied. Note that in Fig. 2-3. that 
while changing the development shifts the location of peak spectral energy deiLsity. it keeps the 
area under the curve (and thus the total energy) constant. 

Ochi Six Parcimeter Spectrum 

The Ochi six parameter spectrum is a generalization of the Bretschneider two parameter spec- 
trum in two respects: 

• two spectral peaks are present, each with their own modal frequency and significant wave 
height, and their contributions add together 

• a shape factor. A, can be applied to each spectral peak to skew the curve right or left, 
independent of the significant wave height or modal frequency. 

In the form given in Eq. 2.8. each parameter is an independent variable. However, by 
analyzing 800 wave spectra from the North Atlantic [2], one can obtain a relationship for the 
six parameters in terms of the measured significant wave height. Ochi reports the results of this 
analysis for a family of spectra having a 95% confidence, however the results given in Table 2.2 
represent a single 'Iriost probable'’ spectrum [2]. The model implemented the most i)robable 
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Figure 2-3: Bretschneider S{cv) for Varying Degrees of Develoj)inent 
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Figure 2-4: Bretschiieider S{u) for Varying Sea Severity. Hs 
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Figure 2-5: Ochi Six Parameter “Most Probable” Spectrum 

form for obtaining the six parameters of the Ochi spectrum. A sample plot of the Ochi most 
prol)al)le spectrum, showing its characteristic biniodal shape, is provided in Figure 2-5. 



2.2 Coordinate System Description 

The coordinate system adopted for the model is a two-dimensional (x-z) system whose origin 
translates in the x direction at the antenna tow speed U . For the case of head seas, U is defined 
to be negative and the waves travel in the positive x direction. Conversely, for the stern sea case 
V is positive and the waves again travel in the positive x direction. Figures 2-6 and 2-7 depict 
the coordinate system employed and the bctsic antenna configuration for both the head and 
stern sea cases. The equations listed in Table 2.1 apply to a stationary coordinate system. 
The traiisformation to a moving coordinate system employs the simple relationship 

x' = X -\-Ui (2.11) 
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Figure 2-6; Coordinate System; Head Seas 



Z 




Figure 2-7; Coordinate System; Stern Seas 
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Figure 2-8: Closeup Showing Key Antenna Position Variables 
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where 



.r' is the earth-fixed (stationary) coordinate system 

.r is the moving coordinate system whose origin translates with velocity U (as shown in 
Figure 2-6 and 2-7). 

2.3 Forces Acting on Antenna 

Th(’ vertical forces acting on the antenna are given in the right hand side of Eq. 2.12 

7 oO ^ 

- T{x,t)-^ 4- Cd(u- - v)\w - ?;| +h{xj) + C M A Pwater ^ 



where 

7/7 is the mass per unit length of the antenna, assumed uniform 

V is antenna vertical velocity (i.e. 

is the tension along the length of the antenna, arising from skin frictional drag 
between the submerged portion of the antenna and the water 

(j is the z coordinate of the antenna central axis 

Co is the normal drag coefficient of a cylindrical antenna section as it moves in the z 
direction through the water 

u' is the vertical fluid velocity 

6(.r, t) is the buoyancy force, which is the resultant of its weight (downward force) and water 
displaced (upward force) 

Cm A is the coefficient of added mass for a cylindrical section in a flow field normal to its 
long axis 

Pwater i^ density of water (for seawater of 15° C, Pwater ~ 1025.9 kg/m^) 

7' is the antenna radius 

El is the product of the antenna cross sectional modulus of elasticity and area moment of 
inertia (thus a measure of lateral stiffness). 

Thus, the first term on the right hand side of Eq. 2.12 represents a restoring force charac- 
teristic of a string displaced from its (taut) equilibrium condition. The tension in the antenna 
is calculated from the time varying wetted surface area along the length of the antenna. A 
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tangential drag coefficient of Ct = 0.0035 was used [ 8 ]. (Reynold’s iiurnbers varied from 32.940 
to 494.130 for speed/diaineter combinations of 3 knots/1 inch to 9 knots/5 inches, respectively.) 
The second term represents the iionnal drag force a cylindrical element would exi)erience as it 
is moved through the fluid with its long axis perpendicular to the direction of motion. A value 
of Co = 1-0 was used (assumes infinite fluid domain). The fourth term represents the added 
mass effect exi)erienced by a cylindrical object as the fluid around is accelerating normal to the* 
long axis of the cylinder. A value of Cm a — 2.0 was selected (cLSSunies infinite fluid domain), 
consistent with standard fluid dynamic texts [9]. The final term represents the restoring force 
due to the stiffness of the material, and its distribution with respect to the neutral axis. At 
this preliminary stage of the antenna design, a material similar to Tygon® tubing is planned 
for the outer anteima housing. At a nominal design of 3 inch outer diameter, 0.25 inch wall 
thickness, and modulus of elasticity E" = 1.5 • 10^ Pa, we have a relatively flexible antenna of 
about El = 300 Pa 

It is recognized that some of these forces have been approximated, namely in that the values 
(juoted al)ove are for Co and Cm a in an infinite fluid domain. The antenna oj)erates in close 
proximity to the free surface, and occasionally portions rise completely above the surface of 
the water. A rigorous treatment of the behavior of Co and Cma near the free surface was 
not attempted. Rather, a straightforward linear approximation was employed which varied the 
coefficients Co and Cma from their infinite fluid domain values when the antenna axis was at 
a dei)th of r (i.e. just submerged), to a value of zero when the antenna axis was at a height of 
r (i.e. just breaking out of the water). This dependence of Co and Cma on relative antenna 
depth y is shown in Eqs. 2.13 and 2.14. 

Coin) = + 9 ) (2-13) 

Cma{v) = + -^) (2-14) 



2.4 Discretization 

The antenna, of length L, is subdivided first into lengths of the physical radio antenna elements. 
Tlu'se antenna elements, designed at Lincoln Laboratory, are called “jelly roll” antennas because 
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Conductor: 
Isometric View 



Antenna Element: 
Section View 





Figure 2-9: '‘Jelly RolF' Antenna Element Construction [12] 

they are a conducting sheet rolled up about 1^ times. (See Figure 2-9 for th(' construction of 
an antenna element.) The element length is designated Leiem> This first level of subdivision is 
not only logical but necessary for the exposure statistics calculated by the model. As it turns 
out. this type of antenna element performs well as long as it has some small amount (e.g. 1 
mm) of its outer housing diameter exposed along its entire length. But if the element becomes 
covered by even a small amount of water at any point along its length, antenna performance 
drops dramatically. Thus it is necessary to monitor numerous points within each element to 
determine whether a given element is exposed or not. This leads to the final level of aiiterma 
subdivision: each element is meshed into some integer number M of sections with length As. 
The subdivision process is shown in Figure 2-10. 

An interesting feature of the program version which calculates the sea elevation and fluid 
velocities by Fast Fourier Transforms (FFTs) is that the frequency meshing (Ao;) of the spectral 
energy density function S{uj), is now related to the spatial meshing (Ax) of the sea, and thus the 
spatial meshing of the antenna (As). This means that there is an extra constraint to consider 
in selecting the finest antenna mesh size As; or conversely, given a mesh size As, constraints 



25 



Antenna 

Length 

L 



Element 
^ Length ** 

^:lem 



Neienf- "7 






Element Length 

Leiem 



> 



^ ^elem^:lem 



Lclem= mas 



Mesh 

Length"* 

As 



M = 6 



Figure 2-10: Antenna Snbdivsion 



26 






are imposed on the Ao.’ meshing of the spectral energy density function S{cj). Some freedom is 
regained if one decouples the length of the antenna and the length of the sea simulated. This 
decoupling is implemented in the model, so that for a typical run the antenna may be 25 in 



long with a mesh size of 0.0625 m. while the length of sea modeled is 512 in long with a mesh 
size of 0.0625 m. The key point is that the mesh length between the antenna (As) and the sea 
(A.r) must match. 

2.5 Computational Method 

2.5.1 Sea Domain Properties 

The sea domain variables of interest are 

• sea elevation, li 

• fluid horizontal velocity, v 

• fluid vertical velocity, lu. 

Two methods were employed in arriving at the sea domain properties, each with certain 
advantages and disadvantages. The first method will be called Sine Wave Superposition (SWS). 
and the second method will be called Fast Fourier Transform (FFT). The source code for both 
of these computer programs can be found in reference [10]. 

Sine Wave Superposition (SWS) Method 

As previously ineiitioned in Eq. 2.1, the sea elevation can be obtained by a superposition of 
sine waves having random phase and amplitude prescribed by Eq. 2.2. In like manner the fluid 
velocities can be calculated by the superposition method. The formulas for /i, u and w in an 
earth-fixed reference frame are provided in Eqs. 2.15 - 2.17. 



N 




(2.15) 



i=\ 

N 




(2.16) 
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(2.17) 



N 

w{x,f) = ujsm{kj{ij) x' — cjjf ^ (pj) 

7 = 1 

By invoking the coordinate transfonnation and the dispersion relation (Eq. 2.11 and Table 2.1, 
respectively), one arrives at the forms appropriate to the coordinate system moving with th(' 
antenna. These are presented in Eqs. 2.18 - 2,20. 



h(x, t) 
u{x. t) 
u’(x, t) 



N 

Uj (cj) cos 

7 = 1 



k,{u) X - {cJr 




)f -f 0, 



N 

7=1 

N 

7 = 1 



c^cos 



cjsin 



k^{iu) X — {lUj 
k2{u) X — {cj-i 



UJ^U 

9 

UjfU 

9 



)t + (pj 

)f + (p, 



(2.1») 

(2.19) 

( 2 . 20 ) 



The spectral energy density function S{cj) used in this formulation was a single sided one. The 
advantages of the SWS method are 



• simplicity in coding fluid property subroutines 

• good resolution for S{uj) at lower frequencies, where most of the wave eiierg}^ resides 

• frecpiency meshing (Acj) is unrelated to the spatial meshing (Ax), allowing both antenna 
element length and wave energy frequency resolution to be arbitrarily chosen. 



Fast Fourier Transform (FFT) Method 

This method takes advantage of the fact that Eqs. 2,15 - 2.17 are a Fourier series, and thus the 
functions /i, u.and w are the inverse Fourier transform of some expression. Since the calculations 
to determine the sea domain properties occur at a fixed point in time, the transform variables 
are not lj and /, rather they become b and x. Table 2.3 summarizes the analogous relationships 
between the time-frequency domains and the space-spatial frequency domains. (The first two 
lines may be taken as definitions; the last line utilizes the dispersion relation.) Because the FFT 
algorithm used for the model assumed the presence of both positive and negative frequency 
components, it was necessary to use a two sided spectrum for this method. In order to retain 
the same amount of energy (i.e. area under the spectral energy density curve) as the transition 
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b= 1/L 


u; = 27t/ 
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II 


k = {2nf)‘^/g 


u = \/27: gb 



Table 2.3: Comparison of Fourier Transform Domain Relationships 



from one-sided to two-sided spectra is made, it is necessary that the formula for the coni})onent 
aiii})litudes be adjusted from that given in Eq. 2.2. In addition, it is convenient to incor])orate 
the random phase o, into a complex component amplitude q. With these changes in mind, the 
formulas j)roviding the sea doniain properties are given in Eqs. 2.21 - 2.25. 



Re exp i{2^xb^x + {27^b^U — \/2Txgbi) /.)j | 

Re c^{bi) yj2ixgbi exp i{2TibiX -h {2'xbiU — \/2T^gbi) t)j 



, 7=1 

N 



, 2=1 



i S{hi) Ab cos((p,) 
Im{c,(6J} = sin{4>,) 

Then the Fourier transform pairs become 



where 



and 



(2.21) 

( 2 . 22 ) 



h{x) 
u{x) 

ic{x) = Re ^ —i Ci{bj) >/27rg6j exp i{2'xbiX + {2-KbiU — \/2Trgbt) f) ^ (2.23) 

Re{ci(6,)} 



(2.24) 

(2.25) 



H{b) <=> k{x) 


(2.2G) 


\/2ii gbi H {b) u{x) 


(2.27) 


—i\/27:gbt H{b) <=> w{x) 


(2.28) 


N 


y^c,(6)exp {2nbtU - \/2ngbij 
2 = 1 


(2.29) 


H{b) ^ h(x) 


(2.30) 
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is defiiieci to mean 



N 

H{b) = ^^Ii{x)exp{—i 27 ibx) (2.31) 

1=1 

1 ^ 

H{b) exp {i 27 ibx) (2.32) 



By setting the complex component amplitudes c^ to obey the relation 



Ci{b^) = o,{-bi) (2.33) 



that is. the positive frequency components are the complex conjugates of the negative frecpiency 
components, then the imaginary parts of Eqs. 2.21-2.23 are zero, and the Re{ } can l)e dropj^ed 
from the equations. When converting between expressing the spectral energy density function 
as a function of uj to another variable, such as b, the basic relations of Eqs. 2.34 and 2.35 are 
useful. 



S{ij) du 
S{b) 



S{b) db 
du 



(2.M) 

(2.35) 



An important feature of the EFT method is the linkage of the spatial frequency domain 
meshing (A5), and the sea spatial meshing {Ax). As previously mentioned, the key to the FFT 
method is to make the antenna spatial meshing. As, match the sea domain spatial ineshing, 
Ax. This is accomplished by using the relations given in Eqs. 2.3G - 2.39, and ensuring that 
the sea Ax equals the desired antenna As. 



Ab = 


1 

j — ^min 

^sea 




(2.3G) 


^max ~ 


^seaAb A^sea 

2 ~ 2Lsea 


1 

~ 2 Ax 


(2.37) 


^min ~ 


\/2tT gbrnm 




(2.38) 


^max ~ 


\/27r^75max 




(2.39) 



Th(^ primary advantage of the FFT method is computational speed. The FFT method 
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Component 


Frequency, u;, 
SWS Method 


FFT Method 


1 


0.347 


0.347 


2 


0.381 


0.491 


3 


0.414 


O.COl 


4 


0.448 


0.G94 


5 


0.482 


0.776 


G 


0.515 


0.850 


7 


0.549 


0.918 


8 


0.583 


0.981 


9 


0.G17 


1.041 


10 


0.650 


1.097 


Notes: 


1. FFT case: 
components 


LU G [0.347,22.2] with 4096 frequency 


2. SWS case: 
components 


uj G [0.347,24.0] with 700 frequency 



Table 2.4: Comparison of Frequency Resolution for the SWS and FFT Methods 

executes approximately 10 times faster than the SWS version. A further increase in execution 
speed is accomplished by not computing the FFT every time step, At. Rather, the sea data 
( /?, u, ic) is computed at nAt intervals, and linearly interpolated in between. The integer n is 
determined by the computer model using the highest sea energy frequency component, 
and the initial estimate of the time step. At. This ensures that the linearly inter])olated sea 
evolves in a smooth and continuous fashion. 

Its disadvantage is that it has poorer low frequency resolution of the sea wave energy than 
its SWS -counterpart. This is demonstrated by a typical run case in Table 2.4 where the lowest 
10 fr(Kiuencies defined are listed for both methods. As one can readily calculate from integrating 
Eq. 2.7 with respect to u. for a sea severity of 2.0m significant wave height, 67.7% of the sea 
wave energy lies below a frequency of 1.097 rad/sec. Thus the FFT method suffers in that it 
])laces most of the frequency components at high frequencies, where little of the wave energy’ 
resides. This undesirable result is a consequence of choosing a b space meshing with constant 
Ab, which is necessary to compute FFT’s in the b and x domains. (Constant Ab meshing 
produces a non-constant Auj meshing with Aui = \/27Tgbin\n — \/f — 1^.) 
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2.5.2 Integration Method 



The time marching integration method implemented is an explicit forward Euler scheiiK?. Such 
schemes are characterized by using information at time t to calculate functions at time t — A/. 
A simple example of how the forward Euler method is used is given below. 

Solve 



dll 



= V 



dt 

n(0) = 1 



(2.40) 

(2.41) 



The tiiiK' derivative is approximated using the Euler method as 



du 



_ 



which for constant time step At becomes 



du 

It 



t=t” 



At 



having adopted the notation of Celia and Gray [11], where 



(2.42) 



(2.43) 



C" = u(t") 



(2.44a) 



and t" is the time step. Then the solution proceeds step by step, beginning with the initial 
value of u given in Eq. 2.41 

= w(0)2 = l (2.45) 



du 

df 



t=o 



then solving Eq. 2.43 for we obtain 



t/i=[/0 + At ^ 

dt 



t=o 



(2.46) 



With the value of u now known at time t (C/^), we can calculate 



du 

~df, 



^=1 ^ ^ 



(2.47) 
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and 



1/2 = + At 




(2.48) 



111 like iiianiier, we move forward in time to calculate any desired 

Spatial derivatives ( alternate between forward centered and backward centered as th(‘ 
derivative order n increases, accomplishing a space centered scheme. Thus the algorithm can 
be classified as a forward time, centered space (FTCS) method. For iiiany designs considered, 
the antenna was approximated as a flexible string of negligible stiffness. With this in mind. Ecj. 
2.12 i)econies a string vibration type equation with damping. As such the Couraiit condition 
[11. p. 231] becomes a useful starting point in selecting a time stej) size. The Courant condition. 
Eq. 2.49 



c At 
A.t 



< 1 



(2.49) 



where 



r is the speed at which waves (i.e., disturbances) propagate down the antenna 

At is the time step size 

Ax is the antenna spatial mesh size 

is used to calculate an initial time increment based on the string wave velocity given by Eq. 



2.50. 




(2.50) 



In Eq. 2.50 

T is the string tension, and 
p is the linear string density. 

Because the tension varies with time, and the maximum tension varies with sea severity 
(due to tangential drag associated with the square of local fluid velocities), the calm water 
tension is used in Eq. 2.50 to form an initial estimate of the time step. Then by multiplying the 
calm water time step by a safety factor, a time step appropriate to the sea severity is arrived 
at. (For a 3.0 m significant wave height, this safety factor is typically set to 0.25) When the 
antenna is modeled as having stiffness, the velocity at which disturbances propagate down the 
antenna is now frequency dependent, and the above implementation of the Courant condition 
no longer applies. Numerical experimentation was used to determine the safety factor which 
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I^roduced an acceptable time step size. 



2.5.3 Boundary Conditions 



The hoiiiidary conditions imposed on the aiitemia endpoints are given in Eqs. 



9(0) = h{0) 
q{L) = h{L) 



(2.51) 

(2.52) 



That is, the antenna follows the sea surface at its eiidi)oints. This is envisioned to be acconi- 
l)lish('d at the trailing end by a drogue having i)Ositive buoyancy, and at the tow i)oint l)y one 
of several options. These options include 

• a lifting body 

• a surface following buoy 

• a leader of positive buoyancy 

• some combination of the above features. 

The problem pertaining to the antenna leader (that portion from the submarine tow point up 
to the surface) is being studied by another graduate student at MIT. and will not be addressed 
in this thesis. 

2.5.4 Initial Conditions 

The antenna simulation begins at t = 0 with 



q{x) = h{x) 

dq{x) dh{x) 

df dt 



(2.53) 



(2.54) 



that is. the anteima is placed at the sea surface, arid given an initial vertical velocity matching 
that of the sea surface. 
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Chapter 3 



Verification of Model Seas 



3.1 Significant Wave Height 



Th(' raiidoiii seaway was constructed by the SWS and FFT methods, as described in Sections 
2.5. land 2.5.1. By running the sea generation subroutines, and sampling the wave amplitude 
record, one can determine the significant wave height of the sample. This can then be compared 
to the significant wave height which was originally specified. From 0chi[2] we have that the 
sample's significant wave height can be found by Eqs. 3.1 and 3.2, 





1 ^ 

= -Yhl 

n - 1 ’ 

1=1 


(3.1) 


H, 


- 4.01 • (7 


(3.2) 



where 

hi is the sampled wave amplitude, measured from the calm water level 
a is the standard deviation of the sampled wave amplitudes 
Hs is the measured significant wave height of the sample. 

(The “ symbol denotes that the quantity is computed from a sanqile, and therefore may 
dcwiate somewhat from the theoretical value.) 

A sea severity of 3.0 in significant wave height was simulated by both the SWS and the 
FFT methods. Twelve wave amplitude records were analyzed by sampling the seas at 8.0 m 
intervals. The results are summarized in Table 3.1. Good agreement to the specified significant 
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k: 



Sample 


Mean(/ij) 

(ill) 


FFT 

a 

(m) 


Hs 

(m) 


Mesin{hj) 

(ill) 


sws 

a 

(m) 


Hs 

(ill) 


1 


-0.0008 


0.7535 


3.0214 


-0.0005 


0.7463 


2.9928 


2 


0.0002 


0.7512 


3.0124 


0.0002 


0.7498 


3.0068 


3 


-0.0003 


0.7490 


3.0037 


-0.0006 


0.7464 


2.9932 


4 


-0.0003 


0.7481 


2.9997 


-0.0006 


0.7494 


3.0051 


5 


0.0006 


0.7511 


3.0119 


-0.0010 


0.7475 


2.9976 


6 


0.0010 


0.7516 


3.0140 


0.0010 


0.7484 


3.0009 


7 


0.0001 


0.7504 


3.0092 


0.0000 


0.7503 


3.0088 


8 


-0.0001 


0.7509 


3.0112 


-0.0002 


0.7501 


3.0081 


9 


-0.0006 


0.7490 


3.0033 


0.0018 


0.7483 


3.0007 


10 


-0.0004 


0.7489 


3.0032 


-0.0011 


0.7525 


3.0170 


11 


0.0002 


0.7492 


3.0044 


0.0002 


0.7479 


2.9991 


12 


0.0002 


0.7500 


3.0075 


-0.0005 


0.7505 


3.0095 


Std. Dev.’s: 


0.0005 


0.0015 


0.0061 


0.0008 


0.0018 


0.0072 


Averages: 


0.0000 


0.7502 


3.0085 


-0.0001 


0.7489 


3.0033 



Tal)le 3.1: Siimilated Significant Wave Heights for a Specified Sea Severity of 3.0 in 



wave lieight is achieved, with an average overshoot in of 0.28% and 0.11% for the FFT and 
S\\"S methods, respectively. 

3.2 Wave Period 



Deteniiiiiatioii of the wave period (or wavelength, frequency, etc.) becoiries difficult for a 
random seaway, where wave shapes and heights are constantly evolving. However, it is quite 
necessary to verify the wave period, because it is possible to correctly model the sea's significant 
wave height, but incorrectly model the wave period (which the research inadvertently proved). 
The method used to verify the wave period invoked the superposition principle: the sea waves 
generated by both methods are ultimately a superposition of individual sine waves, thus the 
verification of wave period was restricted to seas represented by a single harmonic wave. An 
iini)ortant aspect of this process is to verify that the antenna sees the correct “encounter fre- 
(lueiicy," i.e. an effective change of the wave frequency as the antenna is towed in a seaway. 
The encounter frequency, is given by 



ujI-U 

CJe = ^0 

9 



(3.3) 
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Figure 3-1: Variation of Encounter Frequency vs. Tow Speed 

where 

L^’o wave frequency in a stationary coordinate system 

U is the anteniia tow speed 
g is the acceleration due to gravity. 

Thus for towing into head seas { U negative), the encounter frequency goes up, and for 
towing with stern seas ( U is positive), the encounter frequency goes down. The eiicouiiter 
frequency was measured for a harmonic wave of frequency = 1.0 rad/sec at tow speeds from 
-G.O ni/s (towing into head seas) to 6.0 m/s (towing with stern seas). The resulting j)lot of uJe 

ji 

vs. Jj is shown in Figure 3-1. The plot of vs. U is expected to be linear with a slope of — ^ 
(se(^ E(j. 3.3), which for ujo = 1.0 rad/sec has a value of -0.1019 ni“F The measured slope was 
-0.1017 ni“\ which shows the proper dependence of encounter frequency on speed has been 
captunHl. 
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Chapter 4 



Exposure Statistics 



The primary goal of the effort was to calculate exj)osure statistics as various design parameters 
were varied. However, before proceeding on to the exposure statistics results, the issue of 
*‘what constitutes an exposed (i.e., receiving) antenna element?” must be addressed. The 
experimentally observed performance characteristics of the “jelly roll” antenna elements drove 
this definition. It was found that for a level element in calm water the antenna gain is stable 
right up to the point at which the outer skin becomes completely submerged. It was also 
found that if any portion of the outer skin was covered by water (such as might happen in a 
wave field), then the gain of the element suffered drastically. Thus an exposed element was 
defined to be one which was exposed by some small threshold height (set to Q.Q A Diameter for 
all simulations) at every point along the length of the element. In a discrete computer model 
some approximation is necessary with regard to “every point along the length of the element.” 
Thus each antenna element was oversampled by an integer factor, M, which the user prescribes. 
Typical designs considered had an element length of 0.5 m, and were oversampled at M = 8. 

4.1 Selection of Simulation Series 

111 order to bound the problem, it was necessary to discuss with the customer (U.S. Navy) what 
the operational profile recpiired, and with designers (Lincoln Laboratory) what the engineering 
aspects of the antenna gain required. These discussions led to the following conclusions on the 
design sjiace to be simulated: 
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tow speed: U < 10 knots 



• 

• sea severity: Hg < 3.0 in (sea state 4 or less) 

• length: vary about L = 25 m 

• diameter: vary about Diameter = 3 inches 

• density: vary about p = 400 kg/rn^ 

• antenna element length: vary about Le/em = 0.5 m. 

In addition to performing simulation series to support the above design space, it was neces- 
sary to devote some simulations to verify that the results are, in fact, meaningful and consistent . 
Toward this end simulation series were made to 

• compare exposure statistics results of the FFT and SWS methods 

• vary duration of the simulation to determine what constitutes an adequate length of time. 

It is not claimed that these simulation series represent an exhaustive list of all the parameters 
which could be varied, nor is it claimed that the model is completely verified by the consistency 
runs made. However the trends presented will be valuable for the design engineers as they 
consider performance trade-offs, and will allow them to better explore the design space of 
interest. 



4.2 Results 

4.2.1 Statistics Collected by the Model 

The computer model was programmed to track three key statistics: 

• average exposure - the average fraction of antenna elements exposed throughout the entire 
duration of the simulation 

• threshold met - the average fraction of the time that some specified number of antenna 
elements (or more) were exposed 
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Table 4.1; Exposure Data Matrix File Format 



• average count erdetectioii length - the average length of the antenna completely exposed 
(i.e.. completely out of the water). 

(The counterdetection length, while not in the original list of desired statistics, was added 
because of the tactical significance of compromising submarine stealth.) 

In addition, a record is made at each sample time of the condition of every antenna element, 
such that when the simulation is finished, an exposure matrix is available for further ])ost- 
jn-ocessiiig. The structure of the exposure matrix file (an ASCII text file) is shown in Table 4.1. 
(A "1" indicates an exposed antenna element, a “0” indicates a non-exposed element. The .r, 
represent antenna element locations, the t, represent sample times.) Examples of information 
which could be extracted from the exposure matrix include 

• average drop-out time of a specific element 

• average (spatial) length of wave wash-overs 

• histograms and standard deviations of the number of exi>osed elements 

• exposure statistics for an array not fully j)opulated (i.e., one could examine only the 
columns .t, which had antenna elements installed) . 

Three outi)Ut files are generated for each run of the computer model. They are the 

• End-state .summary file - lists input parameters used for the simulation, forces and veloc- 
ities for the final time, and a summary of statistics collected 
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Parameter 


Description 


Method 


SWS or FFT used to generate seas 


Hs 


significant wave height 


Devel 


sea development factor 


^min 


lowest wave frequency component 


^max 


highest wave frequency component 


^ freq 


number of (positive) wave frequency components 


L 


length of antenna 


^elem 


number of antenna elements 


Le/em 


antenna element length 


M 


element oversample factor 


A.S 


antenna mesh size 


Diameter 


antenna diameter 


Pant 


antenna density 


El 


antenna stiffness 


Duration 


simulation duration 


U 


tow s[)eed 


^thresh 


threshold number of elements required for speci- 
fied array gain 



Table 4.2: Simulation Parameter List 



• Exposure data matrix file - provides exposure condition of each element at each specified 
time increment (refer to Table 4.1) 

• Snapshot file - provides a visual display of the antenna as it floats on the sea surface. 
These output files are described in greater detail in Appendix A. 

The input parameters which are si)ecified for each run of the computer model are listed in 
Table 4.2. 

In the- exposure statistics results given below, each series of simulations is presented in the 
following manner: 

• thc^ parameter input table used in the series is listed 

• plots are provided of the three exposure statistics tracked 

— average fraction of elements exj)osed 

— fraction of the time that the threshold number of elements were exposed 

— average counterdetectable length of the antenna 
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It is important to realize that these plots represent averages in time and space. This is 
necessary in order to display the results as functions of key antenna design parameters (length, 
diameter, etc.), which is most useful from a designer’s viewpoint. From the output files, it is also 
possible to analyze how the statistics vary as time progresses, or how the statistics vary down 
th(' array from the tow point to the trailing end. Because of the large number of simulation 
runs, only a few representative plots of this type of information will be shown. 
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Parameter 


Description 


Method 


sws 


Hs 


various 


Devel 


1.0 


^min 


0.4 rad/sec 


^max 


24.0 rad/ sec 


^ freq 


700 


L 


25 111 


^elem 


50 


^elem 


0.5 111 


M 


' 8 


As 


0.0G25 111 


Diaineter 


3.0 inch 


Pant 


O 

00 


El 


0 Pa in"^ 


Duration 


180 sec 


U 


various 


^thresh 
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Table 4.3: Series 1 Input Parameters 



4.2.2 Series 1: Statistics vs. Sea Severity and Tow Speed 
Input 

The input parameters for Series 1: Statistics vs. Sea Severity and Tow Speed, are given in 
Table 4.3. The tow speeds used were 3. 5, 7, and 9 knots; the significant wave heights used 
were 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 meters. 

Results 

Average fraction of elements exposed. Figure 4-1 shows a rather weak dependence of 
the average fraction of elements exposed versus significant wave height, decreasing as the sea 
severity increases. However, the plot shows a strong dependence of average exi)osure versus 
tow speed. 

Fraction of time threshold met. With a required threshold of 40 out of 50 elements. 
Figure 4-2 again shows that tow sj)eed is the dominant factor in the exposure probabilities.. 
The slight rise in threshold met versus significant wave height for the 7 and 9 knot data is an 
interesting feature, which will be discussed further below. 
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Average counterdetection length. As one would expect. Figure 4-3 shows the average 
length of the antenna totally exposed increases with tow speed and sea state. The tow speed 
contributes added tension in the antenna, and the wave height provides the troughs necessary 
for a section to become exposed. 

Histogram Analysis for Number of Exposed Elements 

As mentioned above, there is slight rise in ’threshold met’ for 7 and 9 knot runs as the sea 
s('Verity increases. This is not obvious at first, but can be exi)lained as follows. Figures 4-4 - 
4-6 show the histograms of the number of exposed elements over time running at 7 knots in 
significant wave heights of 1.0, 2.0, and 3.0 meters. For low sea states, the waves have little 
energy, and thus have little ability to affect the exposure statistics. This leads to a relatively 
tightly grouped distribution, a.s in Figure 4-4. Because the mean is below the threshold required, 
the 'threshold met' is very low. (As in Figure 4-2.) But as the significant wave height increases, 
the waves contain more energy, and cause greater variation in the exposure statistics. This 
causes the distributions to flatten out. while the mean exposure stays relatively constant. As 
the upper tail of the distribution moves outward from the mean, it provides more “hits’' above 
the required threshold of 40 elements. This trend is confirmed in Figures 4-5 and 4-6. 

4.2.3 Series 2: Statistics vs. Length 
Input 

The in|)ut parameters for Series 2: Statistics vs. Length, are given in Table 4.4. The threshold 
nuniber of elements was varied in proportion to the length (and the number of elements), but 
was kept at 80% of the total nuniber of elements. Element length was kept constant. 

Results 

Average fraction of elements exposed. The average fraction of elements exposed drojis off 
(piickly at first, then becomes nearly flat, as shown in Figure 4-7. This effect must be studied 
more' closely, since the boundary conditions have a much greater relative inqiact on a short 
antenna than a long one. (Recall the tow point and trailing end are set to float at the sea 
surface.) 
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Average Exposure 



0.950 
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Figure 4-1: Average Exposure vs. Tow Speed and Significant Wave Height 



Parameter 


Description 


Method 


sws 


Hs 


2.5 m 


Devel 


1.0 


^'min 


0.4 rad/sec 


^’max 


24.0 rad/sec 


^ freq 


700 


L 


various 


^elem 


various 


^elem 


0.61 m 


M 


10 


As 


0.061 m 


Diameter 


3.0 inch 


Pant 


408 kg/m^ 


El 


0 Pa in'* 


Duration 


180 sec 


U 


5 knots 


IV 

^thresh 


various {O.SNelem) 



Table 4.4: Series 2 Input Parameters 
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Threshold Met 



1.000 

0.900 

0.800 
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0.000 




0 0.5 1 1.5 2 2.5 3 3.5 



Significant Wave Height (m) 



-©—Tow Speed 3 kts — b— T ow Speed 5 kts —a— T ow Speed 7 kts Tow Speed 9 kts 



Figure 4*2: Threshold Met vs. Tow Speed and Significant Wave Height 
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Average Counterdetection Length (m) 
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Figure 4-3: Counterdetection Length vs. Tow Speed and Significant Wave Height 
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Number of Exposed Elements 



Figure 4-4: Histogram of Exposed Number of Elements at 7 kts; 1.0 m 
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Figure 4-5; Histogram of Exposed Number of Elements at 7 kts: Hs = 2.0 m 
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Figure 4-6: Histogram of Number of Exposed Elements at 7 kts: Hg= 3.0 m 
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Parameter 


Description 


Method 


FFT 


Hs 


2.5 111 


Devel 


1.0 


^min 


0.347 rad/sec 


^max 


22.206 rad/sec 


^ freq 


4096 


L 


25 m 


^elem 


various 


^elem 


various 


M 


various 


As 


0.0625 m 


Diameter 


3.0 inch 


Pant 


408 kg/rn^ 


El 


0 Pa m"^ 


Duration 


180 sec 


U 


5 knots 


^thresh 


various (0.8 Ag/em) 



Table 4.5: Series 3 Input Parameters 



Fraction of time threshold met. Because the threshold was kept at 80% of the total 
iiiiiiiber of elements, we see a sharp drop off in the threshold met versus length in Figure 4-8. 
If the threshold were maintained at some constant number of elements (e.g., 15), one would see 
the opposite trend as length increased. 

Average counterdetection length. The rise in counterdetection length in Figure 4-9 is a 
combination of two factors: maximum tension grows proportional to the length, and as antenna 
length grows beyond the average wavelength it bridges the wave peaks more often. 

4,2,4 Series 3: Statistics vs. Element Length 

Input 

The input parameters for Series 3: Statistics vs. Element Length, are given in Table 4.5. The 
total antenna length was held at 25 meters, and the element lengths used were 0.25, 0.5. 1.0, 
and 5.0 meters. The oversample factor M was varied to keep the antenna mesh length As 
constant at 0.0625 meters. 
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Average Exposure 




Antenna Length (n 



Figure 4-7: Average Exposure vs. Length 
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Figure 4-8: Threshold Met vs. Length 
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Counterdetection Length (m) 




Figure 4-9; Counterdetection Length vs. Length 
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Results 



Average fraction of elements exposed. The results in Figure 4-10 show the effect of 
requiring that a contiguous length (the element length) be exposed. As the contiguously exi)osed 
length reciuireiiient grows, the probability of meeting it decreases. 

Fraction of time threshold met. Figure 4-11 shows that the threshold met dro]>s drasti- 
cally a^ element length grows. The reason can be traced back to Figure 4-10. which shows th(' 
mean exi)osure falls l)elow the threshold of 80%. The more the average exposure falls l)(4ow the 
recpiired threshold number of elements, the more rapidly the ’threshold met’ droi)s off toward 
zero. This pattern would be expected for any distribution with central tendencies. (“Central 
tfuidency'* distributions are distributions which have a higher probability of returning a valiK' 
near the mean than values far away from the mean. The normal (bell-shaped) distribution is a 
classic example of a distril)utioii with a central tendency.) 

Average counterdetection length. Varying the element length should have no effect on 
th(' totally exposed length (the elements are assumed to have flexibility given by El). Figure 
4-12 supports this claim. 

4.2.5 Series 4: Statistics vs. Dicuneter 
Input 

The input i)arameters for Series 4: Statistics vs. Diameter, are given in Table 4.0. All para- 
meters were held constant, except for the antenna outer diameter, which took on values of 1.0, 
2.0. 3.0, and 5.0 inches. 

Results 

Average fraction of elements exposed. The average fraction of elements exposed is a 
strong function of antenna diameter, as seen in Figure 4-13. This is to be expected, since 
antenna tension and normal drag forces (per unit length) grow proportional to the diameter, 
and buoyancy forces (per unit length) grow as diameter squared. 

Fraction of time threshold met. (Figure ??) Typical for distributions with central 
tendencies, and a mean which varies as in Figure 4-13. 
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Figure 4-10: Average Exposure vs. Eleiiient Length 
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Threshold Met 




Figure 4-11: Threshold Met vs. Element Length 
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Counterdetection Length (m 
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II 



Par 2 uneter 


Description 


Method 


FFT 


Hs 


2.5 in 


Devel 


1.0 


^’min 


0.347 rad/sec 


^’max 


22.20G rad/sec 


^ freq 


409C 


L 


25 111 


^elem 


50 


^elem 


0.5 m 


M 


8 


As 


0.0C25 m 


Dimneter 


various 


Pant 


408 kg/m^ 


El 


0 Pa 


Duration 


180 sec 


U 


5 knots 


^thresh 


40 



Table 4.6: Series 4 Input Parameters 



Average counterdetection length. Figure 4-15 shows a strong dependence of counter- 
detection length on diameter. This is again because the buoyancy forces depend on diameter 
scpiared. and thus large diameter designs “ride” the waves, instead of cutting through them. 

4.2.6 Series 5: Statistics vs. Density 
Input 

The input i)ararneters for Series 5: Statistics vs. Density, are given in Table 4.7. All parameters 
were held constant, except for the antenna density, which took on values of 200, 300, 408, 500, 
600, and 700 kg/m^. (The density of the prototype antenna element is 408 kg/m^.) 



Results 

Average fraction of elements exposed. The result shown in Figure 4-16 looks very much 
like an inverse of the result from the previous series. This is because the downward weight 
force (per unit length) grows proportional to the density, while other forces modeled remain 
unaffecte^d as the antenna density varies. 
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Average Exposure 




Figure 4-13: Average Exposure vs. Diameter 



Parameter 


Description 


Method 


FFT 


Hs 


2.5 m 


Devel 


1.0 


^’min 


0.347 rad/sec 


^max 


22.206 rad/sec 


^ freq 


4096 


L 


25 m 


^elem 


50 


L/elem 


0.5 in 


M 


8 


As 


0.0625 m 


Diameter 


3.0 inch 


Pant 


various 


El 


0 Pa in'^ 


Duration 


180 sec 


U 


5 knots 


^thresh 


40 



Table 4.7: Series 5 Input Parameters 
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Figure 4-14: 
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Counterdetection Length (m 




Figure 4-15: Counterdetection Length vs. Diameter 
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Fraction of time threshold met. (Figure 4-17) Typical for distributions with central 
tendencies, and a mean which varies as in Figure 4-lG. 

Average counterdetection length. Figure 4-18 shows a knee in the curve at the 408 
kg/ni‘^ density. While other factors shift this knee left or right, the density of the ju'ototype 
antinma element happens to be about the lowest value one would desire from a counterdetection 
stand])oint at this speed-sea state combination. So from a connectivity viewpoint, lower density 
is desirable; from a stealth viewpoint, lower density is not desirable. Such trade-offs occur 
routinely in submarine operations. 

4.2.7 Series 6: Statistics vs. Duration 

This series and the one following are intended to support the credibility of the above simulation 
results. Th(‘ purpose of this validation series was to demonstrate that the 180 second duration 
used in the computer runs was long enough to be representative of a random process. There 
were two conflicting goals with regard to the duration: 

• computer run time dictated the shortest duration feasible 

• statistical confidence dictated the longest duration possible. 

Thus three runs were conducted at three different durations (for a total of nine runs), in 
order to understand the variance involved in the statistics presented in the first 5 series of 
simulations. 

Input 

The input i)araiiieters for Series 6: Statistics vs. Duration, are given in Table 4.8. Three runs 
each were made for duration values of 60, 180, and 720 seconds. 

Results 

Referring to Figures 4-19 - 4-21, each duration has three bars plotted which correspond to the 
three runs made for that duration. The dashed lines represent the best estimate of the actual 
statistic plotted, since it is the average of the three 720 second run results. Above each duration 



63 




Figure 4-lG: Average Exposure vs. Density 
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Figure 4-17; Threshold Met vs. Density 
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Figure 4-18: Counterdetectioii Length v.s. Density 
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Parameter 


Description 


Method 


FFT 


Hs 


2.5 m 


Devel 


1.0 


^min 


0.347 rad/sec 


^max 


22.20C rad/sec 


^ freq 


409C 


L 


25 ni 


^elem 


50 


^elem 


0.5 m 


M 


■ 8 


As 


0.0C25 m 


Diameter 


3.0 inch 


Pant 


408 kg/m^ 


El 


0 Pa ni^ 


Duration 


various 


U 


5 knots 


^thresh 


40 



Table 4.8: Series 6 Input Parameters 



value is the mean and standard deviation associated with those three simulations. As expected, 
the longer the simulation lasts, the smaller the standard deviation becomes. The goal of the 
computer model is to show trends which will assist the designers in choosing the optimum 
anteniia parameters. At this stage of the simulation effort, the precise values for element 
exposure cannot be predicted without more extensive experimental validation efforts. Thus the 
standard deviation associated with the 180 second duration was judged to be a satisfactory 
compromise between computer run time and statistical confidence. 

4.2.8 Series 7: Statistics vs. Sea Generation Method 
Input 

This series ran identical parameters (within the limits of the meshing conditions) using the 
SWS and the FFT method, in order to ensure the results were consistent between the methods. 
Table 4.9 shows the values used, and indicates which parameters were not matched exactly 
l)ecause of frequency meshing issues. The significant wave height was set to 1.0 and 3.0 meters, 
in order to compare the methods over the full range of sea states simulated. 
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Average Exposure 
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Figure 4-19: Average Exposure vs. Duration 
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Threshold Met 
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Figure 4-20: Threshold Met vs. Duration 
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Counterdetection Length (m 



0.06 



0.04 



0.02 



Mean =-0.088 
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Mean = 0.095 
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Figure 4-21: Counterdetection Length vs. Duration 
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I 




■mt I ^ f4l* "Dfcl ; H 




Parcuneter 


Description 


Method 


FFT / SWS 


Hs 


Various 


Devel 


1.0 


^min 


0.347 / 0.4 rad/sec 


^max 


22.206 / 24.0 rad/sec 


^ freq 


4096 / 700 


L 


25 m 


^ el em 


50 


Lelem 


0.5 rn 


M 


8 




0.0625 ni 


Diameter 


3.0 inch 


Pant. 


CO 

00 

o 


El 


0 Pa m'' 


Duration 


180 sec 


U 


Various 


^thresh 


40 



Table 4.9: Series 7 Input Parameters 



Results 

Figures 4-22 - 4-25 compare the results between the FFT and SWS methods. The agreement 
between the methods is very good; the differences are on the order of the standard deviations 
found in the previous section. 
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Figure 4-22: Coiiiparison of SWS and FFT Methods: Average Exposure for 1.0 m Hg 
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Figure 4-23: Comparison of SWS and FFT Methods: Average Exposure for 3.0 Hg 
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Figure 4-24: Comparison of SWS and FFT Methods: Threshold Met for 1.0 in Hs 
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Figure 4-25: Comparison of SWS and FFT Methods: Threshold Met for 3.0 m 
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Chapter 5 



Future Efforts 



Th(' Buoyant Cable Array Antenna ATD is funded for FY 99 and FY 00. Should this proof- 
of-concei)t ])hase show that the floating antenna principles are feasible and cost effective, the 
Navy is expected to fund the development of a design which can stand up to the rigors of 
fleet deployment. Lincoln Laboratory is scheduled to test a prototype antenna in June 1999 
off the coast of Hawaii. Data will be collected on the electromagnetic performance of the 
antenna elements under varying conditions of source strength, source elevation angle, and source 
azimuthal angle. In addition, the testing will document sea state and towing speeds so that 
perforniance can be related to the sea environment. 

Work at MIT will continue to support the floating antenna effort. Another graduate student, 
under the supervision of Professor Jerome Milgram, will continue and extend the scope of the 
simulation effort. Areas which will be investigated to improve the capability of the computer 
model include 

• 2 dimensional modeling of the sea surface 

• arbitrary angle between tow heading and direction of seas 

• more detailed modeling of the hydrodynamics at the free surface 

• higher order integration schemes 

• validation against at-sea test data. 
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Appendix A 



User’s Guide to Antenna 



The User's Guide is broken down into two sections: 

• Input File 

• Output File. 

Section A.l will go through the Antenna_FFT input file line by line, and include program 
features or limitations associated with each parameter. At the end, features unicpie to the 
Antenna^ SWS input file will be addressed. Section A. 3 will briefly describe the contents of each 
of the three output files, and recommended methods of viewing or extracting the information. 
The output files are the same for both program versions. 

A.l Antenna _ FFT Input File 

A. 1.1 Header 

The input file for Antenna^ FFT is shown in Appendix B. The first 7 lines are for informational 
purposes, and have no effect on the program. However, note that all variables are given in SI 
units, and that there must be 7 lines prior to the first variable defined. 

A. 1.2 Run Identifier {mn^id) 

The run_ id is a 7 character variable used to name the 3 output files. The 3 output files are 
named 
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• sn’// run_id 

• ed’// run _ id 

• 'xl’// Tun_id. 

(The // represents FORTRAN90 string concatenation.) Output files are located in the local 
directory from which the program is run. A sample pattern for run_id is ran_ zrf='001_FFT\ 
Characters which Matlah® would interpret as operators (/, -f , *. etc.) should hv avoided in 

run_id if the Matlab® output o])tioii is selected. 

A. 1-3 FFT Power of 2 {nnf) 

Thv FFT algorithm is designed to operate on arrays which have an integer power of 2 elements. 
nnf sets this power of 2. making 

N.sea = 2^^f (A.l) 



where 

A^-sea is the number of sea frequency components and spatial locations. 

The number of sea modeled in the spatial domain is generally much higher than the antenna 
length requires, however the speed advantage of the FFT method more than outweighs the 
penalty of carrying along the extra sea information. 

A. 1.4 Number of Antenna Elements (num of elem) 

This is the number of physical antenna receiver elements in the array, assuming that 

• the array is fully populated 

• the elements have zero space between them. 

A. 1.5 Antenna Element Length {element length) 

This is the length of each antenna element. 
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A. 1.6 Antenna Oversample Factor (M) 

Each antenna element is broken down into an integer ninnber, M, of sections. These sections 
correspond to the antenna spatial rnesh size As. By enforcing the condition that As = Ax. th(‘ 
computer model constructs the spatial frequency meshing. Ah. such that 



, _ 1 
max 2 ^5 


(A.2) 


L-sexi = As N sea 


(A.3) 


Ah = — 

L.sea 


(A.4) 



By invoking the dispersion relation, we convert from the b domain to the lj doniain. using 

(jj = y/271 gb (A. 5) 

This process shows how the spatial meshing relates to the sea energy frequency doniain 
ineshing. As one can see, the fact that Nsea is constrained to be a power of 2 introduces some 
limitations on the mesh size As, since it is required that 

element Jength = M As (A.G) 

A. 1.7 Maximum Number of Samples {max num of samples) 

This refers to the maximum number of “snapshots” taken for the Matlab® viewer, or the 
exposure data matrix. max_ num_ of_ samples is used to dynamically allocate array storage for 
variables which are sampled periodically in time. 

A. 1.8 Water Density (rho^water) 

Water density in which the antenna operates. 

A. 1.9 Antenna Density {rho^ antenna) 

The average antenna density, which is assumed to be constant over the entire array. This is 
an approximation in that the prototype unit built by Lincoln Laboratory has slightly higher 
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.... r%. .t.v- '. 




density at the antenna element, and slightly lower density between elements. 

A.1.10 Antenna Diameter {diam_ antenna) 

This is the outer diameter of the hose-like structure encasing the antenna. Element performance 
(gain) is controlled by whether this outer diameter is submerged; the diameter of the inner 
conducting slot antenna affects the rotational invariance of the element, and the iiiaxiniinii 
th(H)ietical gain of a fully exposed element. 

A. 1,11 Drogue Size {drog^area) 

Th(' drogiK’ frontal area (in m^) determines the tension at the trailing end of the antenna, 
according to 

T = ^PxvaterU'^ drog.area CD.n (A. 7) 

A value of 0.01 was used for all simulations. 

A. 1.12 Threshold Antenna Exposure {reqd^ant^ exposure) 

As ineiitioned earlier, the antenna gain remains steady until the element just becomes sub- 
merged. In order to capture this performance, and introduce margin for capillary waves not 
modeled, the parameter reqd_ ant_ exposure has been defined. The reqd_ ant_ exposure is th(‘ 
height above the water which must be inaintained by the antenna to be considered exi>osed. 
The reqd_ ant _ exposure is given as a fraction of the antenna diameter. See Figure A-1 for a 
diagram representing an antenna diameter of 1.0 and a reqd_ant_ exposure of 0.04. (0.04 is the 
value used on all simulation results presented.) 

A. 1.13 Tow Speed (C7) 

Th(' tow speed of the antenna is assumed constant with respect to time. Tow speed must be 
greater than zero (in order to generate tension); direction of tow is specified with the next vari- 
able. head_seas. For the case of stern seas only, sea frequency components above omega-crit, 

omega^crit = ~ (A. 8) 
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Figure A-1: Sample Threshold Antenna Exposure of 0.04 

ar(' zeroed out to prevent antenna folding. If omega.crit is less than 1.5 rad/sec. then a warning 
is issued that excessive frequency components are being zeroed out. and program execution 
stops. 



A. 1.14 Direction of Seats Option (head-seas) 

The sea direction is limited to pure head seas or pure stern seas. head_ seas is a logical variable: 
.TRUE, indicates towing into head seas { U < 0), and .FALSE, indicates towing with stern seas 
( U > 0). The sign of U is assigned internally by the program; U should always be positive 
when specified in the input file. 

A. 1.15 Normal Drag Coefficient (CD n) 

This represents the drag coefficient for a cylinder in cross flow. A constant value of 1.0 has 
been assumed. 

A. 1.16 Tangential Drag Coefficient (CD_t) 

This is the skin friction associated with flow along the long axis of the anteniia. A relatively 
smooth outer surface has been assumed in choosing this value, which has been assumed constant. 
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(0.0035 used in all simulations.) 



A. 1.17 Added Mass Coefficient {C^ra) 

The added mass coefficient for a cylinder in cross flow has been set to 2.0. 

A. 1.18 Simulation Time {duration) 

Th(‘ total length of time simulated by the program. While many factors affect the actual tiiiK' 
necessary to perform the siiiiulation (antenna length, number of mesh points. CPU speed, etc.), 
the ratio of actual time to simulation time is roughly 5:1. (For the SWS version, the ratio is 
roughly 50:1.) 

A. 1.1 9 Time Margin {margin) 

Th(' time step size used by the program, t_incr, is given by 



and 

T is the string tension 

is the antenna density 

D is the antenna diameter. 

The tension used initially to estimate the time step is the calm water (local fluid i)article 
velocity, u — 0) tension at the tow point. When dynamic _ tension = .TRUE., then the actual 
tension varies with x and t as the actual wetted surface area changes, and as the local fluid 
l)article velocity, u, evolves along the length of the antenna. A sample plot of maximum tow 
l)oint tension vs. significant wave height using the parameters of Table A.l is j)rovided in Figure 
A>2. As can be seen in the figure, the maximum tow point tension in a 3.0 m significant wave 



tSncr = margin — 



As 



(A.9) 



where c is the string wave propagation velocity 




(A.IO) 
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Maximum Tow Tension (N) 



Parameter 


Description 


Method 


s\vs 


Hs 


various 


Devel 


1.0 


^’min 


0.4 rad/sec 


^’max 


24.0 racl/sec 


N jreq 


700 


L 


25 in 


\T 

^ ' elem 


50 


L^elem 


0.5 m 


M 


8 


An 


0.0625 m 


Diameter 


3.0 inch 


Pant 


408 kg/in'^ 


El 


0 Pa ni^ 


Duration 


180 sec 


U 


5 kts 


^thresh 


40 



Table A.l: Tension Input Parameters 
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Figure A- 2; Variation of Maximum Tow Point Tension vs. Significant Wave Height 
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height sea is about triple that of the calm water value. This would lead to a required value* for 
inm'gin of or about 0,58. In practice, a value of 0.25 to 0.33 was used; however, no detailed 
effort was made to determine the maxiinuiii value allowed for margin. 

For the case of dynamic _ tension = .FALSE., the calm water tow tension is also the* max- 
imum tow tension. (Tension is determined entirely by a relative fluid velocity of U . and the 
wetted area corresponding to the calm water submergence depth.) Thus a value much closer to 
1.0 is possible for margin. 

When the antenna is modeled as having lateral stiffness (i.e. flex _ stiff _coef > 0.0), the 
varial)le margin is used to empirically determine a satisfactory time incremeiit. Time step size is 
cpiitc* sensitive to antenna stiffness, and margin must reduced an order of magnitude for typical 
designs considered (e.g. flex _ stiff _coef =100 Pa m"^). 



A. 1.20 Snapshot Output Interval (pic^time) 

This variable allows the user to control the time interval between snapshots, which are recorded 
to file ’sn’/ /run_id. These snapshots record both the sea and antenna elevation data along the 
(uitirc* length of the antenna, at a resolution of As. A Matlab® viewer script file is included 
in Appendix D, which allows the user to stej) through the snapshots frame by frame, pic_ time 
has been set to 1.0 second for all run results presented. 



A. 1.21 Time of First Snapshot (t^first^pic) 

The time to begin recording snapshots is user defined. This is useful to ignore transient behavior 
at the start of the simulation, or to save only a representative sample at the end of a lengthy 
simulation. 



A. 1.22 Statistics Output Interval (stat time) 

This variable allows the user to control the time interval between exposure data matrix records, 
which are recorded to file ’ed ’/ / run_ id. While stat_ time is independent of ptc_ time, a value of 
1.0 second was also used for this variable for all simulations. Selection of a stat_ time excessively 
small results in subsequent records containing essentially the same information, and contributes 
little to the statistics of the sample. 
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A. 1.23 Threshold Number of Elements {num_ of _elem_ thresh) 

This variable allows the user to specify the lowest number of exposed antenna elements necessary 
to achieve the desired design anteima array gain. At each stat_ time, the number of exposed 
elements is compared to num__of_elem_thresh\ the result (a “1” when the threshold is met. 
“0" otherwise) is written at the end of the summary output file. W / j run_id, along with any 
other variables tracked over the course of the simulation. Currently, the variables trackcxl over 
time are 

• tow point tension {tow _tension{ )) 

• iiumber of elements exposed {is_exposed{ , )) 

• threshold met {thresh_met: “O’’ or “1”, as above) 

• couuterdetection length {length _totally_exposed{ )). 

A. 1.24 Tension Model Option{dynarnic_tension) 

(See discussion under Section A. 1.19.) 

A. 1.25 Snapshot Output Oipt\on{matlah_output) 

The siiai)shot output file {''six' // run __id) is optional, since it does not contain any unique sta- 
tistical informatioii. and requires a relatively large amount of disk space. (Average file size is 
about 1.4 MB for a 180 second simulation.) The primary benefit of examining the snapshot 
output file is to learn the physical behavior of the antenna in a seaway. In addition, the output 
is useful to detect when the time step size is too small, which is characterized by ra[)id un- 
I)hysical antenna oscillations. matlah_output=.’TWJYj. means that snapshot outj)ut is enabled, 
.FALSE, means snapshot output is disabled. 

A. 1.26 Significant Wave Height {sig^wave^ht) 

This variable sets the sea’s significant wave height, a i)arameter used to describe the sea sever- 
ity in both the Bretschneider and Ochi sj^ectral models. A})})endix E correlates sea state to 
significant wave height for a fully develoi)ed sea. 
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A. 1.27 Development {development) 

As discussed in Section 2.1.2, this variable allows the user to shift the sea spectral energy modal 
frequency uj) (for developing seats) or down (for decaying seas). This parameter api)lies only to 
th(' Bretschiieider spectral model. 

A. 1.28 Spectral Model { spectral _ model) 

Two spectral model subroutines are available: the two parameter Bretschiieider spectrum, and 
th(' Ochi ‘’niost probable” six parameter spectrum. (‘T” means Bretschiieider, “2” means Ochi.) 
Th(» Bretschiieider spectral model was used for all simulations presented. 

A. 1.29 Antenna Lateral Stiffness {flex_ stiff _coef) 

This variable represents the resistance of the antenna to bending, and is given by 

flex-stiff-coef = El (A. 11) 



where 

is the modulus of elasticity 
1 is the area irioinent of inertia. 

A. 2 Input File Features Unique to Antenna_SWS 

The differences in the Antenna_SWS input file pertain to the discretization of the sea circular 
frequency components. 

A. 2.1 Number of Sea Frequency Components {num of_freqs) 

This variable is the number of discrete sea frequency components in the spectrum. num_ of_freqs 
= 700 was used for the simulations, however as few as about 100 comporients provide adecpiate 
sea representation. The frequency cornpoiients are evenly spaced in the to domain. 
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A. 2. 2 Frequency Coverage {omega^min, omega^raax) 

The freciiieiicy coinpoiients are evenly spaced from 07iiega_min to omega^max. For the 
simulations presented oniega_min = 0.4 rad/sec and omega_max = 24.0 rad/sec. 

A. 2. 3 Spectral Model {spectral _model) 

One extra option is added: for spectral _model = 3. the sea surface is described by a pure 
harmonic wave of amplitude sig_wave_ht and circular frequency omega_max. 

A. 3 Output Files 

A. 3.1 End-State Summary File (’xl’//run_ id) 

This file provides 

• a printout of forces, velocities, and positions at the end of the run 

• a record of the input parameters used for the simulation 

• averaged statistics results, and 

• a time history of variables tracked by the program. 

This file is suitable for viewing with a spreadsheet application, such as Microsoft Excel®. 

A. 3. 2 Exposure Data Matrix File (’ed’//mxn_id) 

This file, whose data structure is given in Table 4.1, contains the raw exposure data information 
for each antenna element at each time sampled. This allows one to analyze statistics across 
time or space. This file is easily analyzed by a spreadsheet application. 

A. 3. 3 Snapshot File (’sn’//mn_ id) 

The snapshot file can be viewed with a user written code, or by Snaps, nr a simple Mat lab® 
routine ]:>rovided in Appendix D. It provides views of the antenna as it floats on the surface at 
two scales: one view is large enough to see the most extreme waves, one zooms in on the sea 
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surface by using the local water elevation as the zero for the z axis. The latter view also shows 
th(' antenna central axis, the antenna upper limit, and the antenna lower limit. 



88 



Appendix B 



Input File for Antenna FFT 



Input File for Antenna_FFT Simulation Program 

SI units throughout; useful conversions follow 
1 inch = 0.0254 m 
1 foot = 0.3048 m 



1 knot = 


0.5144 m/s 




OOO.FFT 


run_id 


7 characters long 


13 


nnf 


# of components. N_sea=2**nnf 


25 


num_of _elem 


# of actual slotted antenna elements 


0.5 


element .length 


physical slotted antenna length (m) . 


8 


M 


# of mesh points/element .length. 


300 


max_num_ of .samples 


to dimension time arrays 


1025.9 


rho.water 


seawater at 15 deg C (kg/m‘3) 


408.0 


rho. antenna 


density of antenna (kg/m“3) 


0.0762 


diam.antenna 


antenna outer diameter (m) 


0,01 


drog.area 


drogue area normal to flow (m‘2) 


0.04 


r e qd_ ant _ expo sure 


threshold vertical exposure for reception, (x/diameter) 


2.5722 


U 


tow speed (use POSITIVE #) (m/s) 


.TRUE. 


head.seas 


TRUE: antenna is towed into head seas 


1.0 


CD_n 


normal flow drag coefficient 


0.0035 


CD_t 


fluid tangential drag coefficient (p. 50 Berteaux) 


2.0 


C_m 


added mass coefficient 


4.0 


duration 


total simulated run time (sec) 


0.27 


margin 


saiftey factor applied to time step size 


1.0 


pic.time 


snapshot interval for Matlab viewer (sec) 


0.0 


t_f irst.pic 


time of first snapshot (sec) 


1.0 


stat.time 


statistics interval (sec) 
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..M 



■s? 



I 



20 


num_of_elem_ thresh 


# of exposed elements reqd for specified gain 


.TRUE. 


dynami extension 


FALSE: steady-state calm water model to find tension 


.TRUE. 


matlab.output 


TRUE: Matlab viewer output is ON 


2.5 


sig_wave_ht 


significant wave height (m) 


1.0 


development 


modal frequency shift factor 


1 


spectral_model 


l=Bretschneider ; 2=0chi 


0.0 


f lex_stif f _coef 


E*I (Pa m"4) 
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Appendix C 



Input File for Antenna SWS 



Input File for Antenna_SWS Simulation Program 



SI units 


throughout; useful conversions follow 


1 inch = 


0.0254 m 




1 foot = 


0.3048 m 




1 knot = 


0.5144 m/s 




OOO.SWS 


run_id 


7 characters long 


25 


num_of _elem 


# of actual slotted antenna elements 


0.5 


element_length 


physical slotted antenna length (m) . 


8 


M 


# of mesh point s/element .length. 


100 


num_of _freqs 


# of discrete freq components 


300 


max_num_of _samples 


to dimension time arrays 


1025.9 


rho_water 


seawater at 15 deg C (kg/m"3) 


408.0 


rho_antenna 


density of antenna (kg/m“3) 


0.0762 


diam_ antenna 


antenna outer diameter (m) 


0.01 


drog_area 


drogue area normal to flow (m"2) 


0.04 


reqd_ant_exposure 


threshold vertical exposure for reception, (x/diameter) 


2.5722 


U 


tow speed (use POSITIVE #) (m/s) 


.TRUE. 


head_seas 


TRUE: antenna is towed into head seas 


1.0 


CD_n 


normal flow drag coefficient 


0.0035 


CD_t 


fluid tangential drag coefficient (p. 50 Berteaux) 


2.0 


C_m 


added mass coefficient 


4.0 


duration 


total simulated run time (sec) 


0.27 


margin 


saftey factor applied to time step size 


1.0 


pic_time 


snapshot interval for Mat lab viewer (sec) 


0.0 


t_f irst_pic 


time of first snapshot (sec) 


1.0 


stat_time 


statistics interval (sec) 
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20 


num_of _elem_ thresh 


# of exposed elements reqd for specified gain 


.TRUE. 


dynamic ^tension 


FALSE: steady-state calm water model to find tension 


.TRUE. 


matlab_output 


TRUE: Matlab viewer output is ON 


0.40 


omega^rain 


lower cutoff freq of sea spectrum (rad/sec) 


15.0 


omega_max 


upper cutoff freq of sea spectrum (rad/sec) 


2.5 


sig_wave_ht 


significant wave height (ra) 


1.0 


development 


modal frequency shift factor 


1 


spectral.model 


l=Bretschneider ; 2=0chi; 3=pure sine wave 


0.0 


flex_stif f _coef 


£♦1 (Pa ra'‘4) 
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Appendix D 



Source Code Listing for 



Snaps.m 



’ snaps displays antenna and sea surface as a series of snapshots 

! displays series of antenna snapshots 

! N=number of antenna elements 
' M=number of snapshots taken 



number of rows in file = M+1 
number of cols in file = 2*N+1 



File format (ASCII) 
col 1 : 

zero (a dummy), time 1 time M 

col 2: 

x(l), z(l) at time 1, ... , z(l) at time M 
col N+1; 

x(N), z(N) at time 1, ... , z(N) at time M 
col N+2: 



x(l), h(l) at time 1, ... , h(N) at time M 
col 2*N+1: 

x(N), h(N) at time 1 h(N) at time M 



clear 

load snOOO_SWS -ascii 



[M ,N] =size (snOOO.SWS) ; 
N=(N-l)/2; 
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t=sn000_SWS(2;M+l,l) ; 
x=snOOO_SWS(l,2:N+l) ; 
z=sn000_SWS(2:M+l ,2:N+1) ; 
h=sn000_SWS(2:M+l ,N+2:2+N+l) ; 
diam=snOOO_SWS (1,1); 
y=z-h; 
for j=l:N 
zz( j)=0.0; 

threshl ( j)=0. 04*diara; 
thresh2 ( j ) =-0 . 04*diam ; 

end 

' find m2Lx/niin z for consistent plot scales 
2 _min=min(h) ; ! z_rain is a row vector 

z_rain=min( 2 _min) ; ! z_min is now a scalar 

z_raax=raax(h) ; 
z_raax=max(z_raax) ; 
f igure(l) ; 
for j=l:M 
If j>l 
pause ; 

end 

subplot (2,1,1); 

plot(x,z(j , ;) ,x,h(j , : ) , x,zz); 

title ( [ 'Time' , blanks (3) ,num2str(t( j))] ) ; 

ylabeK'z, h (m) O ; 

axis([0.0 x(N) z_rain z_max] ) ; 

subplot (2,1,2); 

plot(x,y(j , :) , ’k: ' ,x, threshl , ’b-’ ,x,thresh2, 'b-’ ,x,zz, ’g- 
axis([0.0 x(N) -.05 .05]); 
xlabeK’x (meters) O ; 
ylabeK'y = z - h (ra) O ; 



end 
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Appendix E 



Sea State Data for Fully Arisen Seas 
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Significant 
Wave Height 
(in) 


Modal 

Frequency 

(rad/sec) 


Wind 

Speed 

(knots) 


Sea 

State 


0.1 


3.96 


4.2 




0.2 


2.80 


5.9 




0.3 


2.29 


7.3 


1 


0.4 


1.98 


8.4 




0.5 


1.77 


9.4 




0.6 


1.62 


10.3 




0.7 


1.50 


11.1 




0.8 


1.40 


11.9 


2 


0.9 


1.32 


12.6 




1.0 


1.25 


13.3 




1.1 


1.20 


13.9 




1.2 


1.14 


14.6 




1.3 


1.10 


15.2 


3 


1.4 


1.06 


15.7 




1.5 


1.02 


16.3 




1.6 


0.99 


16.8 




1.7 


0.96 


17.3 




1.8 


0.93 


17.8 




1.9 


0.91 


18.3 




2.0 


0.89 


18.8 


4 


2.2 


0.85 


19.7 




2.4 


0.81 


20.6 




2.6 


0.78 


21.4 




2.8 


0.75 


22.2 




3.0 


0.72 


23.0 


5 


3.2 


0.70 


23.8 




3.4 


0.68 


24.5 




3.6 


0.66 


25.2 


6 


3.8 


0.64 


25.9 




4.0 


0.63 


26.6 




4.5 


0.59 


28.2 




5.0 


0.56 


29.7 




5.5 


0.53 


31.2 




6.0 


0.51 


32.6 




6.5 


0.49 


33.9 


7 


7.0 


0.47 


35.2 




7.5 


0.46 


36.4 




8.0 


0.44 


37.6 




8.5 


0.43 


38.8 




9.0 


0.42 


39.9 





Table E.l: Sea State Data for a Fully Arisen Sea 
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